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Abstract 

We study a simple Standard Model (SM) extension, which includes three families of right-handed neutrinos with generic 
non-trivial flavor structure and an economic implementation of the invisible axion idea. We hnd that in some regions of the 
parameter space this model accounts for all experimentally conhrmed pieces of evidence for physics beyond the SM: it explains 
neutrino masses (via the type-I see-saw mechanism), dark matter, baryon asymmetry (through leptogenesis), solve the strong 
CP problem and has a stable electroweak vacuum. The last property may allow us to identify the Higgs field with the inflaton. 
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1. Introduction 

Although no unambiguous signal of physics beyond the 
SM (BSM) has appeared so far at the LHC, there is no 
doubt that the SM has to be extended. Neutrino oscillations, 
which lead to the existence of small (left-handed) neutrino 
masses, and the observational evidence for dark matter (DM) 
is enough to state that the SM is incomplete. 

Other unsatisfactory features of the SM are an insufficient 
baryon asymmetry of the universe, the strong CP, gauge hi¬ 
erarchy and cosmological constant problems. 

Moreover, precision calculations [1, 2] indicate that the 
SM potential develops an instability at a scale of the order 
of 10'° GeV, for central measured values of the SM param¬ 
eters. This is not particularly worrisome per se because the 
probability of tunneling to the absolute minimum, where life 
is impossible, is spectacularly small [2]. However, it may 
lead to some issues during the exponential expansion of the 
early universe (inflation) [3, 4, 5]. Moreover, the (absolute) 
stability up to the Planck scale Mpi may lead to the possibil¬ 
ity of Higgs inflation [6, 7, 8, 9], linking particle physics and 
cosmology: this is interesting because it provides us with 
relations between particle physics and cosmological observ¬ 
ables. The presence of such an instability in the SM is not 
firmly confirmed because of non-negligible uncertainties on 
the top mass and the QCD gauge coupling; but, if confirmed, 
it would suggests that right-handed neutrinos (at scales suit¬ 
able for the see-saw mechanism and thermal leptogenesis) 
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and the physics of the QCD axion may be relevant for the 
issue of the electroweak (EW) vacuum instability and there¬ 
fore inflation. 

The aim of this paper is to identify a simple and well- 
motivated model where the following signals of BSM physics 
can all be addressed and which adds to the SM only right- 
handed neutrinos and the extra fields needed to implement 
the axion idea: 

1. Small neutrino masses. We adopt the perhaps simplest 
explanation: the type-I see-saw mechanism based on 
right-handed neutrinos. The addition of right-handed 
neutrinos also symmetrize the field content of the SM 
giving to each SM left-handed particle a right-handed 
counterpart. 

2. Dark matter. As a DM candidate we consider the 
axion [10], a light spin-0 particle whose existence is 
implied by the spontaneous symmetry breaking of a 
U(l) symmetry, the Peccei-Quinn (PQ) symmetry [11] 
that explains why strong interactions do not violate CP. 
In particular, we consider the invisible axion model 
proposed by Kim, Shifman, Vainshtein and Zakharov 
(KSVZ) [12], which has a simple structure and a small 
number of free parameters. 

3. Baryon asymmetry. In order to explain such asymme¬ 
try we make use of (thermal) leptogenesis [13], which is 
implemented with the same right-handed neutrinos that 
allow the light neutrinos to have masses. 
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4. Inflation and vacuum instability. As we stated be¬ 
fore, the inflaton could be identified with the Higgs bo¬ 
son provided that the EW vacuum is stable^ taking into 
account energies up to the Planck scale. We therefore 
look for regions of the parameter space where the EW 
vacuum is stable, even for central values of the SM ob¬ 
servables. 

5. Strong CP problem. The solution we consider is the 
first and most famous one; the PQ symmetry, the same 
symmetry leading to the axion DM candidate above. 

It is important to note that the hrst two points represent a 
proof of BSM physics, while the others are indications, al¬ 
though very plausible ones. The spirit here is similar to the 
one of [15], focusing on problems 1, 2 and 3 and adding 
to the SM field content only right-handed neutrinos with 
masses below the EW scale. In this case, indeed, the right- 
handed neutrinos can signihcantly contribute to dark mat¬ 
ter [16] (see also [17] for a review and further references) 
and neutrino oscillations provide a mechanism to generate 
baryon asymmetry through a different version of leptogene- 
sis [18]. Moreover, it is possible to extend this framework 
to include the axion idea and to look for simultaneous solu¬ 
tions of problems 1, 2, 3, 4 and 5. Here there is no claim that 
the simple model we study is the only one able to address all 
these issues. 

In the list above we did not include the gauge hierarchy 
and the cosmological constant problems because they can 
both be addressed with anthropic arguments^ [19]. On the 
other hand, there seems to be no anthropic solution to the 
strong CP problem; thus technical naturalness appears to be 
the only possible way to explain the small value of the QCD 
9 angle. 

Let us summarize now the contents of the article. In sec¬ 
tion 2 we define the model. In section 3 we discuss the ob¬ 
servational constraints on its parameters. The theoretical in¬ 
gredients for the extrapolation up to Mpi are provided in sec¬ 
tion 4. In section 5 we investigate whether the model can 
have a stable EW vacuum taking into account energies up 
to Mpi. One of the conditions for stability is that the Higgs 
quartic coupling remains always positive. There is, however, 
another condition to be fulfilled to ensure a stable vacuum. 


’ By adding non-renormalizable operators with independent coefficients 
one may enter the region of metastability [14], we do not consider this pos¬ 
sibility in the present paper. 

^The gauge hierarchy problem can of course be solved in a technically 
natural way (e.g. with SUSY, composite Higgs, etc) in models that explain 
some of the issues mentioned above [20]; also, SUSY large extra dimensions 
[21 ] olfers a possible way to address the cosmological constant problem; but 
this is done at the price of introducing many more fields than those of the 
model studied here and sometimes the necessity of an ultraviolet completion 
at much smaller energies. 


This section contains the central new results of this paper. 
Einally, in section 6 we provide our conclusions. 

2. The model 

We consider the model with Lagrangian: 

— -^gravity "h .£sM "h -^N "h .£axion. ( 1 ) 

where repeated indices understand a summation. The gauge 
group of the model is the SM one; 

Gsm = SU(3)eXSU(2)iXU(l)y. 

-Cgravity are the terms in the Lagrangian, which include the 
pure gravitational part and the possible non-minimal cou¬ 
pling between gravity and the other fields. In particular, the 
term proportional to \H\^'R, where H is the Higgs doublet and 
7? is the Ricci scalar, plays an important role in Higgs infla¬ 
tion [ 6 ]. XsM is the SM Lagrangian (minimally coupled to 
gravity). X.n is the part of the Lagrangian that depends on 
the right-handed neutrinos A, (i=l,2,3); 

JIn = iNimi + + YijUHNj + h.c.^ , (2) 

where M,y and Ey are the elements of the Majorana mass ma¬ 
trix M and the neutrino Yukawa coupling matrix Y, respec¬ 
tively. Thanks to the complex Autonne-Takagi factorization, 
we take M real and diagonal without loss of generality; 

M - diag(Mi,M 2 ,M 3 ), 

where the M, (i — 1,2,3) are mass parameters, the Majorana 
masses of the three right-handed neutrinos. 

Einally, -Caxion represents the additional terms in the La¬ 
grangian due to the chosen axion model. As stated in the 
introduction, we consider the first invisible axion model (the 
KSVZ modeP). The fields of this model that are not con¬ 
tained in the SM are the following. 

• An extra Dirac fermion. (In Weyl notation) it is a pair 
of two-component fermions qi and q 2 in the following 
representation of Gsm 

qx ~ (3.l)o, qi ~ (3.l)o. (3) 

Namely they form a colored Dirac fermion with no in¬ 
teractions with the gauge fields of SU(2)/, x U(1 )f. 

• An extra complex scalar. This scalar A is charged un¬ 
der U(l)pQ and neutral under Gsm- 


^For a recent interesting work where another axion model and scalar 
generations of neutrino masses are considered, see [22]. 
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The Lagrangian of this axion model is 

2 

-£axion “ ^ E qjlj>qj + \d^A\^ - {yqiAqx + h.c.) - AV(H,A) 
j=i 

and the classical potential of the fall model is 

V{H,A)^AHm^-v^f + AV(H,A), (4) 

where 

AV(H,A) = d^(|A|2 - /2)2 + - v2)(|A|2 - /2). 

The parameters v, fa and y can be taken real and positive 
without loss of generality. The PQ symmetry acts on qi, q 2 
and A as follows 

qi ^ e'^^^qu qi^e'“'^q 2 , A^e~'°A, (5) 

which forbids an explicit mass term Mqqiq 2 + h.c.. The SM 
fields and the right-handed neutrinos are instead neutral un¬ 
der U(l)pQ. Moreover, there is the accidental symmetry 

q\ ^-qi, q 2 ^q 2 , A ^-A. ( 6 ) 

This model has the advantage of being simple and having 
(in addition to the SM and type-1 see-saw parameters) only 
three real parameters: Aha, 4a and y; it is the most general 
one given the field content and symmetries described above. 
In particular, notice that Aha is the only tree-level coupling 
between the axion and SM sectors. 

The EW symmetry breaking is triggered by the vacuum 
expectation value (VEV) v ^ 174GeV of the neutral compo¬ 
nent Hq of the Higgs doublet. After that the neutrinos acquire 
a Dirac mass matrix 

niD - vY, (7) 

which can be parameterized as 

niD = ( niDi , mD2 , '«d3 ) , ( 8 ) 

where moi (i = 1,2,3) are column vectors. Integrating out 
the heavy neutrinos Nj, one then obtains the following light 
neutrino Majorana mass matrix 

mD2m]j2 '«D3'«D3 

my — -1-1-. (y) 

Ml M 2 M 3 

By means of a unitary (Autonne-Takagi) redefinition of the 
left-handed SM neutrinos we can diagonalize my to obtain 
the mass eigenvalues m \, m 2 and TO 3 (the left-handed neutrino 
Majorana masses). Calling Uy the unitary matrix that imple¬ 
ments such transformation, also known as the Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix, that is UfmyUy - 
diag(mi, m 2 , m 3 ), we can parameterize Uy = VyP 12 , where 

/ C 12 C 13 inCo \ 

Vy — j -il2C23 - Cl2il3‘S234'* C12C23 - Si2SiiS23e‘^ Cl3'y23 1 , 

V ■Sl2‘y23 “ ^12^13023^"’ —C 12^23 — S 12^13^23^^ C 13 C 23 / 


with Sij = sm(0ij), Cij = cos( 0 y); O/j are the neutrino mixing 
angles and P 12 is a diagonal matrix that contains two extra 
phases, in addition to the one, 6, contained in Vy: 

/ 0 0 \ 

P12 = 0 0 . ( 10 ) 

Vo 01 / 

Even in the most general case of three right-handed neutri¬ 
nos, it is possible to express Y in terms of low-energy observ¬ 
ables, the heavy masses Mi, M 2 and M 3 and extra parameters 
[23]: 

V 

where 

Dy/Wi = diag( VmT, Vm2, V^), 

= diag(VA^, Va^, y/^) 

and /? is a generic complex orthogonal matrix, which con¬ 
tains the extra parameters. This is useful for us because the 
observational constraints are not directly on Y, but they are 
rather on the low-energy quantities m,, Uy and on M; (see 
section 3). One can show that the simpler and realistic case 
of two right-handed neutrinos [24] below Mpi can be recov¬ 
ered by setting mi -0 and 

/ 0 0 1 \ 

R - cosz - sinz 0 , 

V ^sinz ^cosz 0 J 

where z is a complex parameter and ^ = + 1 . 

The PQ symmetry is broken both spontaneously and by 
anomalies. The spontaneous symmetry breaking is induced 
by fa = (A), leading to the following Dirac mass of {qi,q 2 ): 

M, = yfa. 

Moreover, A contains a (classically) massless particle, the 
axion, which acquires a small mass thanks to the quantum 
breaking of the PQ symmetry, and a massive particle with 
squared mass 

Ml^f]:{AAA+o(^j^y (12) 

As we will review below, the observational bounds imply 
that the corrections O {v^lfa) are very small and will be ne¬ 
glected in the following. 

3. Observational constraints 

We now discuss the observational constraints, which we 
will take into account in the rest of the paper. 
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As far as the neutrino masses m,- (i — 1,2,3) are con¬ 
cerned, data from atmospheric and solar neutrinos tell us 
respectively [25] (see also [26, 27, 28] for previous deter¬ 
minations) 

= 7.50!jS;|^ X 10^^eV^ 

Amh = 2.457!j];j]^:5 X lO^^eV^, 

where Amfj = mf-mj and Am^ = j for normal ordering 
and Aniy = Am ^2 for inverted ordering. 

As far as the mixing angles and phases of the PMNS ma¬ 
trix are concerned, the most recent central values and cor¬ 
responding uncertainties can also be found in [25]; for any 
ordering of the neutrino masses the 3cr ranges are 

0.270 < ^ 0-344, 0.385 < < 0.644, 

0.0188 < si 3 < 0.0251, (13) 

while 6 spans the whole range from 0 to In at her level (for 
example for normal ordering we have while, 

for inverted ordering, d/° = 254^^1). Currently no signiheant 
constraints are known fory 6 i and^ 62 . 

We now turn to the requirements to have successful lep- 
togenesis [13]; neutrinos should be lighter than 0.15 eV and 
the lightest right-handed neutrino Majorana mass M; has to 
fulfill [29] 

M, > 1.7x lO’GeV. (14) 

In order to be conservative we have reported the weak¬ 
est bound, but depending on the assumptions one can have 
stronger conditions"^. Notice, however, that the mechanisms 
of [15] and [18] discussed in the introduction can evade these 
bounds and use right-handed neutrino masses below the EW 
scale; as we will see, it is less challenging to achieve vac¬ 
uum stability in this case. In other models, if the Higgs field 
acquires a large VEV during inflation, [30] argued that the 
subsequent Higgs relaxation to the EW vacuum can generate 
the baryon asymmetry. 

Regarding the axion sector, in order to account for DM 
through the misalignment mechanism [31] (with an order 
one initial misalignment angle) and to elude axion detection 
one obtains respectively an upper (see e.g. [32]) and lower 
bound (see e.g. [33]) on the order of magnitude of the scale 
of PQ symmetry breaking fa'. 

10^ GeV <fa< 10 '^ GeV. (15) 

The upper bound is obtained by requiring that the axion field 
takes a value of order fa at early times, which is what we ex¬ 
pect, but is not necessarily the case; also the precise value of 


^For example if the initial abundance of right-handed neutrinos at T 3> 
Ml is zero then the bound is M/ > 2.4 X 10^ GeV [29]. 


the lower bound is model dependent. Therefore (15) should 
not be taken as sharp bounds, but it certainly gives a plausi¬ 
ble range of /„. Another source of uncertainty is introduced 
if one instead considers light right-handed neutrinos [16, 17], 
which can then contribute to dark matter, as mentioned in the 
introduction; in this case, indeed, the upper bound becomes 
stronger as it is obtained by requiring the axion contribution 
not to exceed the observed dark matter abundance. In any 
case, (15) ensures that /„ » v and the terms O (v^lfa) in 
(12) can be neglected. Moreover, notice that bounds on fa 
can only constrain the ratio Ma! as it is clear from ( 12 ). 
When Ma » v and Mq » v (which we assume) the EW 
constraints are fulfilled. 

In addition to contributing to dark matter, the axion also 
unavoidably manifest itself as dark radiation as it is also ther¬ 
mally produced [34, 35, 36]. This population of hot axions 
contributes to the effective number of relativistic species, but 
the size of this contribution is currently well within the ob¬ 
servational bounds [36]. 

Einally, of course we also have constraints on the SM pa¬ 
rameters. After the discovery of the Higgs boson at the LHC 
[37, 38] the last SM parameter, the Higgs mass, has been 
determined within small uncertainties, and there are no free 
SM parameters anymore. We take the values and uncertain¬ 
ties of the SM masses and couplings given in [2] (see also the 
references therein). The determinations of [2] are not signif¬ 
icantly affected by the presence of the extra heavy degrees of 
freedom. 


4. RGE analysis and thresholds 


Since we want to study the predictions of this model at 
energies much above the EW scale, up to the Planck scale, 
we need the complete set of renormalization group equa¬ 
tions (RGBs). We adopt the MS renormalization scheme 
to define the renormalized couplings and the corresponding 
RGBs. Moreover, for a generic renormalized coupling g we 
write the RGBs as 


dr 




(16) 


where djdr = jf'dldjf and f is the MS renormalization 
energy scale. The y 6 -functions fg can also be expanded in 
loops as 




(47r)^ (47r)"^ 


(17) 


where KAttP" is the n-loop contribution. 

Let us start from energies much above Ma, Mg and Mij. 
In this case the 1-loop RGBs are (see [39, 40, 41, 42] for 
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Figure 1: Phase diagram of the model, showing the region with ab¬ 
solute stability up to the Planck scale. The region where condition 
II is fulfilled is inside the region with Au(jT) > 0. We set the central 
values of the SM parameters at the EW scale and the low-energy 
neutrino parameters; however, we checked that variations of Am^ 
and O^j (within Scr around their central values) and variations of 6, 
Pi and p 2 have a negligible effect on this plot. Moreover, we set 
the lightest neutrino mass m, = 0 , M 2 = lO'^^GeV, M 3 > Mpi and 
z = 0. Switching the sign of f does not change the plot. The axion 
decay constant is set to /„ = iO**GeV and T^fM^) = 0.05. 


previous determinations of some terms in these RGBs) 
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yW - ^gi), 

where g^, g2 and = -\/5j3gY are the gauge couplings 


of SU(3)c, SU(2)/, and U(l)y respectively and y, is the top 
Yukawa coupling. The explicit form of the complete set of 
the RGBs above was not explicitly presented before, but the 
RGBs for a generic quantum field theory (without gravity) 
were computed up to 2-loop order in [43] (see also [44] for a 
computer implementation of them). 

Next, we consider what happens in going from energies 
above Ma to energies below Ma' as discussed in [45,41] one 
has to take into account a scalar threshold effect: in the low 
energy effective field theory below Ma one has the effective 
Higgs quartic coupling 


A^Ah- 



( 18 ) 


This is the result of integrating out the massive scalar degree 
of freedom at tree-level. The reason why this shift occurs 
is because setting the heavy scalar to zero is not a consis¬ 
tent truncation, namely it is not consistent with the equations 
of motion. In practice one should do the following; below 
Ma the RGBs are the ones given above with and 
removed and Ah replaced by A. Above Ma one should in¬ 
clude p^^^ and p^^ and find Ah using the full RGBs and the 
boundary condition in (18) at /i = Ma- 

As far as the new fermions are concerned, following 
[46] we adopt the approximation in which the new Yukawa 
couplings run only above the corresponding mass thresh¬ 
olds; this is implemented technically by substituting Yjj 
YjjOifi - Mj) and y —> y0(p - Mq) on the right-hand side of 
the RGBs. The situation is different from the scalar one, as 
setting the fermion fields to zero below their mass threshold 
is consistent. 

Binally notice that, the SM parameters can run in an en¬ 
ergy range bigger than the one of Y, Aa, Aha and y. There¬ 
fore, we include for them the 2-loop RGB contribution; we 
do not, however, show explicitly the 2-loop part because of 
its complexity. 


5. Stability analysis 

Since we use the 1-loop RGBs of the non-SM parame¬ 
ters, we approximate the Coleman-Weinberg [47] effective 
potential of the model with its RG-improved tree-level po¬ 
tential; we substitute the bare couplings in the classical po¬ 
tential with the corresponding running ones. 

The conditions that ensure the absolute stability of the vac¬ 
uum (Hq) - V and {A) - fa have been studied in [41]: they 
are 
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I. AnipL) > 0 and AAijY) > 0 

II. Ac = AAh{p)Aa{]j) - Ajj^ip) > 0 
















Figure 3: RG evolution of the quartic couplings Ah, Aa and the com¬ 
bination of quartics defined in condition II for stability. The 
vertical solid line indicates the position of the scalar threshold, Ma- 
The stripes on the right indicate the region presumably dominated 
by Planck physics. The values of the parameters are the same used 
in figure 1 . 


Figure 2: The same as in fig. 1, but with a different value of y. 

Notice that once Ah > 0 and A,, > 0 are fulfilled then /l ,4 > 0 
is fulhlled too. The fact that the MS couplings are gauge 
invariant, as proved in [48, 2], guarantees that our results 
will not be affected by any gauge dependence. 

The first condition Ah > 0, at the level of approximation 
we are using, may lead to the possibility of Higgs inflation 
[7, 8, 9]. Therefore having absolute stability may also al¬ 
low us to identify the inflaton with the Higgs field. However, 
one should keep in mind that perturbative unitarity^ is vio¬ 
lated above some high energy scale [50, 51]. Once the back¬ 
ground helds are taken into account, however, the authors of 
[52] hnd that such energy is parametrically higher than all 
relevant scales during the history of the Universe. Never¬ 
theless some extra assumptions on the underlying ultraviolet 
completion are necessary [51, 52, 8]. 

The question of the stability of the EW vacuum has been 
addressed previously in other economic extensions of the 
SM. The SM extended only by adding a single right-handed 
neutrino or three right-handed neutrinos with degenerate 
masses was studied in [46, 40]. Extensions with a singlet 
scalar were considered in [53, 41, 54] and others with one 
right-handed neutrino and an extra real scalar were studied 
in [55]. However, we do not know of any previous work that 
accounted for all problems listed in the introduction®. 


^This unitarity problem can be solved by adding an extra real scalar field 
[49, 41]. The extension of the present analysis to include such scalar is 
beyond the scope of this paper. 

^After posting this article on the arXiv our attention was drawn to the 


In hg. 1 and 2 we show regions of the parameter space 
where the stability conditions are fulhlled for all values of jA 
up to Mpi and others where they are not. The values of the pa¬ 
rameters used in that plot can also explain neutrino masses, 
dark matter, baryon asymmetry and the strong CP problem 
(through the mechanisms discussed in the introduction), ful- 
hlling all bounds of section 3. Moreover, the regions where 
Ah > 0 all the way up to Mpi correspond to the possibility of 
Higgs inhation. In hg. 2 we see that increasing y{Mq) shrinks 
the region where condition II for stability is fulhlled; this is 
because y contributes positively (negatively) to the ranning 
of Aha (Ta), which then increases (decreases) and this makes 
it more difficult to satisfy that condition. We also observed 
that changing the value of Aa(Ma) and fa changes the loca¬ 
tion of that region, so that the size of the parameter space 
that is compatible with absolute stability is larger. Notice 
that hgs 1 and 2 also indicate that lighter right-handed neu¬ 
trino masses favor the stability conditions. This can be qual¬ 
itatively understood: smaller M; generically correspond to 
smaller T,y, Eq. (11), and to a reduced destabilizing effect in 

conditions I and II because of the way Y appears in and 

Ml) 

P^HA- 

In figs. 3 and 4 we show the evolution of the quartic 


interesting Ref. [56]. The authors discuss a model very similar to ours 
and anticipate that all those problems (with the exception of the origin of 
inflation) can be solved: in that work the PQ symmetry is an extension of the 
SM lepton number. This allows to relate the scales fa and M, [57]. However, 
an explicit analysis was not presented in [56]. As we will see now, such an 
analysis here leads to regions where the simultaneous solutions occur and 
others where they do not. 
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Figure 4: The same plot as in fig. 3, but with the values of the pa¬ 
rameters used in figure 2. 

coupling combinations relevant for the stability analysis as 
a function of the renormalization scale. The parameters are 
chosen in a way compatible with the regions of, respectively, 
figs. 1 and 2, where all stability conditions are fulfilled. 
There are no Landau poles below the Planck scale and the 
couplings remain perturbative when the stability conditions 
are fulfilled. The region with stripes on the right corresponds 
to the regime where Planck physics is expected to be dom¬ 
inant; the behavior of the curves there is thus presumably 
unreliable. 

At the same time, it is important to notice that there are 
also regions of the parameter space, where the results on 
the stability analysis obtained in the SM are not significantly 
changed by the addition of A,, qj and A. In the limit Aff^ ^ 0 
the axion sector is decoupled from the rest, and, if the neu¬ 
trino Yukawa couplings are small enough, one recovers the 
SM results at a very good level of accuracy. 

6. Conclusions 

In this paper we have found regions of the parameter 
space of a simple but well-motivated model that can ac¬ 
count for all experimentally confirmed signals of physics be¬ 
yond the SM: neutrino oscillations (through the addition of 
three right-handed neutrinos), dark matter (due to the axion), 
baryon asymmetry (generated by thermal leptogenesis), in¬ 
flation (which could be driven by the Higgs field since the 
EW vacuum can be an absolute minimum for energies up to 
the Planck scale) and the strong CP problem that is automat¬ 
ically solved by the PQ symmetry leading to the axion. 

This model is an extension of the SM, which only adds 
to the SM three right-handed neutrinos as well the scalar 


field and extra colored fermion of the simple invisible axion 
model proposed by KSVZ. 

We have found that there are values of the parameters such 
that the important features listed above are all present to¬ 
gether with perturbativity (always up to the Planck scale). 

An important extension for the present work may be the 
inclusion of quantum gravity, which has been completely ne¬ 
glected here. Some steps in this direction have been taken in 
[42]. But the role of gravitational quantum effects in the sta¬ 
bility issue of the SM is still unclear. 
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